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We construct Heisenberg anti-ferromagnetic models in arbitrary dimensions that have isotropic 
valence bond crystals (VBC) as their exact ground states. The d = 2 model is the Shastry-Sutherland 
model. In the 3-d case we show that it is possible to have a lattice structure, analogous to that of 
SrCu2{BOz)2, where the stronger bonds are associated with shorter bond lengths. A dimer mean 
field theory becomes exact at d — > oo and a systematic 1/d expansion can be developed about it. 
We study the Neel-VBC transition at large d and find that the transition is first order in even but 
second order in odd dimensions. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.50.Ee 



Recently there has been a renewed interest in the 2-d 
Shastry-Sutherland model (SSM) B, owing to its physi- 
cal realisation in SrCu2{BOz)2 [| |uj . The model has an 
exactly solvable ground state. There exists some gener- 
alized antiferromagnetic Hamiltonians with exact ground 
states §, HI. 

SSM was initially constructed as a 2-d generalisation 
of the 1-d Majumdar- Ghosh model 0. Both models have 
a valence bond crystal (VBC) as the exact ground state. 
Other such models, including a 3 dimensional one [|, 
have since been constructed. 

All the abovementioned models can be thought of as 
special cases of the class of models which we will define 
below and refer to as generalised Majumdar-Ghosh mod- 
els (GMGM). Consider Hamiltonians of the form, 

H = JnK (1) 

n 

Where the sum is over all possible triangles formed by the 
sites of the lattice, J„'s are arbitrary positive semidefmite 
couplings and h n are given by, 

h n = SfcJ.Sfo) + S(r i ).S(r Jk ) + S(r fc ).S( ri ) (2) 

S(rj), S(rj) and S(rfc) are the spins at sites r,, Vj and rj. 
respectively and n labels the triangle formed by them. A 
particular can be a part of more than one h n . 

It was noted in Q that if, in the triangle corresponding 
to h n , two of the spins are forming a singlet, then the 
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FIG. 2: Shastry-Sutherland Model 



state will be a ground state of h n . It was also pointed out 
that h n could be more general than given above. It could 
be spin anisotropic, the three terms could have different 
coefficients and the spin could be arbitrary. In suitable 
parameter ranges the dimer state will remain the ground 
state ||. Thus, if it is possible to cover the lattice 
with dimcrs in such a way that each of the triangles that 
appear in H with a non-zero coupling contains a dimer 
then the state with singlets on all the dimers will be a 
simultaneous ground state of all h n and hence that of H . 

The problem of constructing GMGMs then reduces to 
the purely geometric one of assigning the non-zero cou- 
plings to the triangles such that a dimer covering of the 
above type is possible. We will now give a class of solu- 
tions to this problem which naturally generalises the SSM 
to arbitrary dimensions. For simplicity, we will work with 
h n as in equation (|) and with s = 1/2. 

We will first construct the 3-d model and then gener- 
alise to arbitrary dimensions. We first set up the VBC 
and then build the Hamiltonian around it. We take a 
simple cubic lattice and choose the dimers to lie along 
the body diagonals. The body diagonals are assigned as 
follows. The sites are denoted by, x = XL=i Xpi&n- x n' s 
take integer values and e M are three orthogonal unit vec- 
tors. The spin S(x) is paired to the spin S(y(x)) where, 



FIG. 1: A generalised Majumdar-Ghosh Chain 



y(x) =x + D(x) 



(3) 



FIG. 3: The 3-d lattice containing 18 sites, depicting the 
diagonals along which dimers are formed 
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D(x) = £(-1)**+^ (4) 

where we define X3+1 = x\. Two such body diagonals 
are shown in FIG. g. 

Note that D(y) = — D(x) as it should be, since if S(x) 
is paired with S(y), then S(y) should be paired with 
S(x). All the four body diagonal directions occur in equal 
numbers and the VBC has cubic rotational symmetry. It 
is not parity invariant, the other parity choice being given 
by replacing (— l) a; f+ 1 by (— l) 2 ^- 1 in equation (Q). 

We now choose the triangles with nonzero couplings 
as follows. Equation (Q) uniquely associates a body di- 
agonal and hence a unit cube with every dimer, (x, y). 
We give non-zero couplings to the six triangles formed by 
these two sites and each of the other six sites that belong 
to the cube as illustrated in FIG. |I| Thus every such 
triangle has one edge, one face diagonal and the body di- 
agonal containing the dimer. This construction ensures 
that every triangle has one dimer (along the body diago- 
nal) and hence is the exact ground state of Hamiltonians 
of the form given in equation (|l]). The Hamiltonian can 
be explicitly written as, 

3 

where, 

ft(x,/i) = S(x).S(z(x, M ))+S(z(x, M )).S(y) 

+ S(y).S(x) (6) 
z(x, M ) = x+t-rf^ (7) 

z(x, n) and x form the 3 edges emanating from x in the 
direction of the body diagonal D(x) and y is given by 



FIG. 4: A cube containing diagonal bond. All the bonds are 
shown 

equation^). J(x, fi) is the coupling associated with the 
triangle formed by x, y(x) and z(x, ji). 

Consider the simplest case when all the couplings are 
equal, i.e. J(x, //) = J. The triangle corresponding to 
/i(x, jj) contributes a strength J to the edge it contains. 
Each edge is contained in exactly one triangle. Thus all 
edges have bond strengths J. Each triangle contributes 
a strength J to the face diagonal that it contains. Half 
the face diagonals are contained in exactly one triangle 
each and the other half are not contained in any triangle. 
Thus half the face diagonals have bond strength J and 
the others have no bonds. FIG. |] illustrates the situation. 
Finally, each triangle contributes J to the body diagonal. 
Half the body diagonals are contained in six triangles 
each and hence have bond strengths 6J and the other 
half have no bonds. See FIG. [| 

The generalisation to higher dimensions, d = 4, 5, ... is 
straightforward. Simply replace 3 by d in all formulas 
from equation to equation (Q). All the 2 d_1 body 
diagonal directions occur in equal numbers in the d di- 
mensional VBC. The model is a simple hyper-cubic lat- 
tice with bonds of strength J along all the edges and 
along one of the face diagonals of each (d — 1) dimen- 
sional face. There are also bonds of strength 2dJ along 
one of the body diagonals of half the hyper-cubes. The 
construction ensures that the VBC is the exact ground 
state of the model. 

It can be seen that the model reduces to SSM at d — 2, 
shown in FIG. [| The diagonals are given by equation(||). 
The strength of the bonds along the diagonals is 4 J. d — 
2 is a special case in that the (d — 1) dimensional face 
diagonals are also the edges. Thus the strength along the 
edges is 2 J, ie. the strength of edge bonds are half that 
of the diagonal ones. Thus we have recovered the 2-d 
SSM. 

Now we come back to the d = 3 case. As it stands, it 
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FIG. 5: A deformed cube. Now the strongest bond is between 
nearest neighbour sites 

is not very physically realistic since the stronger bonds 
are between spins further apart. The same is true in the 
case of the 2-d model. However, the structure of mag- 
netic ions in SrCu2{BOz)2 can be obtained from the 
original theoretical lattice by moving the sites along the 
diagonals that have bonds. The squares containing the 
diagonals deform to rhombi and the body diagonal con- 
taining the dimer becomes shorter than the edges. As 
we will see, the procedure generalises to 3-d and we can 
obtain an analogous structure where the stronger bonds 
have shorter bond lengths. 

We move all sites along the body diagonals with non- 
zero bond strengths. The new sites are then, 

R(x)=x+^D(x) (8) 

where s is a parameter and D(x) is given by equation(||). 
The cube shown in FIG. ^ deforms to the rhombohedron 
shown in FIG. Before the deformation the sites formed 
a simple cubic lattice. The Hamiltonian was however 
only symmetric under translations by two units. After 
deformation, the lattice periodicity is also halved. It re- 
mains a cubic lattice but with 8 sites in a unit cell. 

The lengths of the edges, face diagonals and body diag- 
onals can be computed and are plotted as functions of s 
in FIG. ||. As we can see, for s more than around .7 or so, 
the edges and face diagonal with bonds become almost 
equal in length, are longer than the body diagonal with 
the bond and shorter than the other face diagonals and 
body diagonals. When s=l the body diagonal becomes 
of zero length. The rhombohedron is then squashed into 
a hexagon lying in the plane orthogonal to the body di- 
agonal. 

We will now examine the models away from the ex- 
act ground state point at d — > oo. We put the bond 
strengths along the body diagonals equal to cZJd/2, along 
edges equal to Je/2 and along the face diagonals to 



FIG. 6: The various bond lengths as functions of s. l e is the 
length of the edges, Ifi that of the face diagonals containing 
bond and I /2 corresponds to face diagonals without bond. Idi 
is the length of the diagonals along which dimers are formed 
and Idi, ld3 and Idi correspond to body diagonals without 
bonds. 

Jf/2. Actually, the VBC is an exact eigenstate when 
J E = J F = J' for all Jd 0- It is proven to be the 
ground state when J' < 0.5 Jd- However, the ground 
state is Neel ordered for Jp = = Jd ■ So there will be a 
phase transistion somewhere. The location and nature of 
the transition at d = 2 have been topics of much activity 

©Hi! linn. 

Since the VBC is the exact ground state in some pa- 
rameter range for all d, it is clear that mean field theory 
in terms of the spin variables fails even as d — > oo. This is 
because one of the bond strengths grows as d and so the 
interactions cannot be approximated by an average field. 
However, if we take the dimers, i.e the two-spin systems 
on the body diagonals to be the basic units, then each 
dimer interacts with ~ d of the dimers around it with 
bonds of strength ~ 1. Thus the mean field theory in 
terms of the dimer variables is exact at d — > oo in this 
class of models. Perturbation theory around the mean 
field Hamiltonian then yields a systematic 1/d expansion 
for the fluctuations. We will now use this mean field 
theory to explore the physics at large d. 

We label the spins as S/ Q , where / labels the posi- 
tions of the centers of the dimers and a = 1, 2 label the 
two spins that form the dimer. The Hamiltonian can be 
written as, 

H = ^ — (S/l + Sj 2 ) 2 + - Jla.jpSla-Sjfj (9) 

/ la, J0 

Where Ji a ,jp denote the edge and face diagonal bonds. 
The mean field Hamiltonian is, 

Hmf = ^2 — ^— (S/i + S/ 2 ) 2 + ^2 Jiajpbia-S.jp 

I Ia,Jf3 

- 2 XI J i*,jpb Ia .b Jl3 (10) 
The self consistency equations are then, 

b/ Q = {Sic) (11) 
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In different parameter regimes, the VBC as well as a 
variety of other phases are possible. The 1/d expansion, 
which is valid in all the phases can be used to analyse the 
phase diagram. In this paper we study the Neel- VBC 
transition at large d. There is a qualitative difference 
between odd and even dimensions and we treat them 
separately. 

In even dimensions, the dimer lattice is bipartite and 
both the sites of a dimer have the same parity. The parity 
of a dimer is defined to be same as that of its sites. The 
Neel-state is described by the ansatz, b/ Q = Pibz, where 
Pj is +1 on one sublattice and -1 on the other. Then the 
mean field Hamiltonian is, 
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(12) 



Where J' = (Je + Jf)/2. For all values of a = J'/Jd, 
all the dimers forming singlets, ie VBC, is a mean- field 
solution with b = 0. The state has energy 0. When 
a > 1/2, the Neel state is a solution with 6=1/2 and 
has energy, 



E Neel = NdJn {1 _ a) 



(13) 



Where N is the number of sites. However, this solution 
has lower energy than the VBC only when a > 1. Thus 
we get a first order transition at a — 1 at d = oo. 

We have also computed the leading order correction 
to the ground state energy in the Neel phase by treating 
(H — Hmf) as a perturbation. We get, 
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At d = 2, the transition now occurs at a s=s 0.8. More so- 
phisticated calculations at d = 2 [jl0| [l^] put this number 
at 0.69. There are indications that the transition may be 
second order p4| or that there may be an intermediate 
phase @ [Of. 

In odd dimensions the two sites of a dimer are not of 
the same parity. By convention, we assign a = 1 for the 
odd site and a = 2 for the even site. Then the Neel 
ansatz is b/ Q = (— l) a bz. The mean field Hamiltonian is 
then given by, 
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(15) 



Where A J = (J E — Jf)/2. In the ground state the dimer 
wave function is given by, 
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(16) 



Where in \l,m), I is the total spin and m the z compo- 
nent and sin9 = 2AJ\h\/( y /J'f ) + (8AJ|b|) 2 ). The self 
consistency equation (p4|) reduces to, 
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(17) 



Thus the transition occurs at Jd = 4A J. The sub-lattice 
magnetization, |b| goes continuously to zero at the tran- 
sition as ~ (1 — ( Jd/4A J)) 1 / 2 . The interesting thing is 
that unlike in even dimensions, the transition point de- 
pends on the difference of Je and Jf- So the VBC can 
occur at relatively low values of Jf>- 

The different physics in the even and odd dimensions 
arises from the fact that the two spins on a dimer belong 
to the same sublattice in former and on opposite ones 
in the latter. Consequently, the dimer wavefunction in 
the VBC and in the Neel state have the same value of 
S z (= 0) in odd dimensions whereas in even dimension 
S z = in the VBC but S z = ±1 (on odd and even 
sublattices) in the Neel state. Since the S z symmetry is 
unbroken in both the phases, the mean field Hamiltonian 
always conserves it. Thus in even dimensions the VBC 
state cannot smoothly transit to the Neel state and we 
get a first order transition whereas in odd dimensions it 
can and we get a second order transition. 

A remark about the scaling of the diagonal bond is in 
order here. For the dimer mean-field theory to be valid 
in both the phases, the diagonal bond has to scale as d, 
so that d scales out of Hmf of equations ( ^2| , |i~5| ). But 
numerically, the critical value of the diagonal bond could 
be small. For example, in odd dimensions the critical 
value of diagonal bond is proportional to A J = ( Je ~ Jf ) 
and can be made arbitrarily small by suitable choice of 
Je and Jf- 

The mean field equations can be solved in the presence 
of an external magnetic field B. The mean- field Hamil- 
tonian in the presence of B = Bz is given by, 



Hmf = — — - (S/i + S/ 2 ) 2 + ^ Jia,j/3bia-Sj@ 

I Ia,J@ 

- BY,S z Ia -\ Ji^jph Ia .hjp (18) 

la Ia,J0 

We concentrate on the regime where the VBC is the 
ground state, in the absence of B. A state in which one 
of the dimers is excited to a triplet with S z — 1 and the 
rest are all in singlets is a self-consistent solution of the 
mean- field Hamiltonian of equation(fl8|), with hj a = 
for those dimers in singlets and hj a = for the one 
in triplet. The energy of this state is dJo/2 — B. Thus 
as we increase the strength of the magnetic field, when 
B > dJu/2, it is energetically favourable to excite as 
many dimers as possible into triplets pointing along B, 
but no two of them being connected by a bond. 

In even dimensions, using the fact that the lattice 
formed by the dimers is bipartite, we shall now show 
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that the maximum fraction of dimers that can be excited 
without any two of them being connected by a bond is 
1/2. 

Let us start with isolated dimers with no bonds con- 
necting them. Now we add edge bonds such that ev- 
ery dimer is connected to one and only one other dimcr. 
This can be done as follows. As mentioned before, both 
the sites of a dimer have the same parity. Each coor- 
dinate of the two sites will be differing by +1 or — 1. 
Choose an even dimer and pick out the site with odd 
xi-coordinate. Put the edge bond from this site in the 
positive xi-direction. This way every even dimer can be 
uniquely connected to one and only one odd dimer. Sup- 
pose these were the only bonds present. Then the con- 
figuration that maximizes the number of triplets without 
two of them being connected is where in every pair of 
connected dimers one is put in a singlet and the other 
in a triplet. Now the original model can be obtained by 
adding the remaining bonds. But in the presence of ad- 
ditional bonds also, each of the pairs that were initially 
connected can at the most have one triplet. Thus the 
maximum fraction of dimers that can be excited without 
any two of them being connected by a bond is 1/2. This 
upper bound can be satisfied by putting triplets on one 
of the sub-lattices and singlets on the other. This state 
is a mean-field solution with b/ Q = for those dimers in 
singlets and hj a = \z for those in triplets. 

The state with the remaining singlets also being ex- 
cited into triplets is also a mean-field solution and has 
energy ^{dJo/2 + dJ' — B) above the half-magnetized 
state. Thus when B > dJu/2 + dJ' all the remaining 
dimers are excited into triplets and the system becomes 
fully magnetized. 

When d is odd, the dimer lattice is not bipartite. We 
still assign parity to dimers as follows. In every dimer one 
site will be odd and the other even. The parity of the 
dimer is defined as the parity of the xi-coordinate of the 
even site. It can be seen that out of the Ad dimers that 
are connected to a particular dimer only 4 are of the same 
parity. Since the critical B, at which it is energetically 
favourable to have isolated singlets, scales as d, this dif- 
ference between odd and even dimensions is insignificant 
at large d. 

For both the cases, there is a plateau at 1/2 magneti- 
zation for dJr> /2 < B < dJo/2 + dJ'. At stronger fields, 
the system is fully magnetized. The 1/2 plateau corre- 
sponds to all the dimers on even sites in the triplet and 
the ones on the odd sites in the singlet state. At d = oo, 
the dimers have only nearest neighbour interactions. Pre- 
vious work |nj [H], [l6| indicates that the other fractions 
are due to longer range interactions induced by fluctu- 



models, d > 2, that are natural generalisations of the 
Shastry-Sutherland model. These are Heisenberg antifer- 
romagnets on hypercubic lattices with bonds of strength 
dJo/2 along half the body diagonals, Je/2 along all the 
edges and Jf/2 along half the (d— l)-dimensional face di- 
agonals. They have isotropic VBCs as their exact ground 
states at Jo > 2Je — Up. A dimer mean field theory 
is exact in the d — * oo limit and a 1/d expansion can be 
developed about it. In this limit, the Neel-VBC transi- 
tion is first order in even dimensions and occurs at a = 1 . 
When leading order (1/d) corrections are included, the 
transition shifts to a = 0.8 at d = 2. In odd dimensions 
the transition is second order and occurs at Jd = 4AJ'. 
In presence of an external magnetic field, B, in both even 
and odd dimensions, the system has a 1/2 plateau from 
B = dJo/2 to dJn/2+dJ' and becomes fully magnetised 
at higher fields. 



At d = 3, we have shown that it is possible to have a 
crystal structure where the stronger bonds have shorter 
bond lengths. Moreover, the mean field theory indicates 
that the transition to the VBC phase occurs at relatively 
small values of Jd ■ Thus a physical realisation of such a 
system seems feasible. 



atiois. Thus it is reasonable that only the 1/2 plateau 



[1 
[2. 
[3 

K 

[5 
[6 

[-: 

[8 
[9 

[10 

[11 

[12 

[13 
[14 

M 



occurs at d = oo. 
To summarise, 
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